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quality engineered for simplicity and 



Plug-in components may be replaced easily and 
quickly for expansion or servicing. 

Non-linear components fit most non-linear row po- 
sitions. Number of configurations is limited only 
by the number of components kept on hand. 
Basic Computer is pre-wired and can be expanded 
simply by plugging in desired components — no 
additional wiring necessary. 

Draws less power than 60 watt bulb. Operates from 
115V, 60 cycle outlet. 

Bus bar power distribution eliminates complex 
cabling and simplifies maintenance. 



Solid state design — instant warm-up — no coolin 
problems. 



Human engineered control panel is inclined for 
easy, finger-tip control. 

Push button potentiometer readout system speeds 
set up — reduces errors. 

Built-in null voltmeter provides direct reading or 
precision null reading. 

Color coded patching modules and accessories pro- 
mote programming efficiency. 



Interchangeable, plug-in components add flexibility, 

make expansion easy 




COEFFICIENT 
SETTING 

POTENTIOMETERS 

for inserting 
equation coefficients 
or problem 
parameters into 
computer. 




NETWORKS 

enable operational 
amplifiers 
to perform 
integration. 




FUNCTION 
SWITCHES 

provide for 

manually 

interchanging 

components 

without 

reprogramming or 
repatching. 




COMPARATORS 

compare a variable 
input voltage to an 
arbitrary bias 
voltage and cause a 
switching operation 
to be performed. 




FUNCTION 
GENERATORS 

electronically 
generate analytic, 
as well as 
arbitrary 
functions of one 
variable. 




REFERENCE 
PANEL 

makes available 
accurate reference 
voltages required 
for problem 
solution. 




TIE POINT 
PANELS 

provide additional 
patch panel 
terminations for 
components inputs 
or outputs. 




MULTIPLIERS 

electronically 
multiply two 
variables of 
either sign. 




OPERATIONAL 
AMPLIFIERS 

are high-gain, 
low-drift, 
chopper-stabilized 
DC amplifiers used 
for addition, 
subtraction, 
integration and 
inversion. With 
other components, 
they also perform a 
variety of 
non-linear 
operations. 



ACCESSORIES 




• Service Shelf (Type 51.039) 

Facilitates maintenance of any 
plug-in computing component 
under normal operating conditions. 




• Feedback Resistor 

Type 646.010 10,000 ohms, ±0.1% wire-round resistor; 
supplied in a BLUE molded plug, is designed for patching 
between summing junction and output terminations of any 
operational amplifier. 

• Input Resistors 

Type 646.006 10,000 ohms, ± 0.1%., an ORANGE banded 
epoxy-encapsulated, wire-wound resistor — male end plugs 
into amplifier summing junction termination, female end 
accepts patch cord plug. 

Type 646.007 100,000 ohms, ± 0.1%, a YELLOW banded, 
expoxy-encapsulated, wire-wound resistor — male end plugs 
into amplifier summing junction termination, female end 
accepts patch cord plug. 

• Resistor Set (Type 5.134) 

Includes the following quantities of Input and Feedback 
Resistors described above: 

15 each Type 646.006 (see above) 

10 each Type 646.007 (see above) 

10 each Type 646.010 (see above) 

• Diode Unit (Type 614.051) 

A WHITE banded, epoxy-encapsulated silicon diode for 
limiting the output of computing components or generating 
non-linear effects. 

• Patchcord Set (Type 5.133) , includes the following: 

10 each Type 510.033-0, 6" long, Color Black 
15 each Type 510.033-1, 12" long, Color - Brown 
10 each Type 510.033-2, 18" long, Color - Orange 
5 each Type 510.033-4, 30" long, Color - Blue 

• Multiple Block (Type 542.605) 

Provides a six hole, off-the-patch-panel tie point for 
interconnecting patch cords or increasing the number of 



Patch Cord output holes of various computing components. 
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BASIC PROGRAMMING PROCEDURES FOR 
TRANSISTORIZED ANALOG COMPUTERS 

INTRODUCTION 

Many of the problems encountered in scientific or engineering endeavors require the solution of mathe- 
matical equations or sets of equations which in most cases are difficult and in other instances are, for 
all practical purposes, impossible to obtain by the classical approach to equation solution. The PACE 
TR-10 Analog Computer aids in overcoming this difficulty by providing the technical worker with a 
simple-to-use computer which permits the rapid solution of linear or non-linear mathematical equations. 

Although the analog machine is correctly termed a computer, it does not perform its computations by 
numerical calculations as does the desk calculator or the digital computer. The analog computer per- 
forms mathematical operations on continuous variables instead of counting with digits. Thus the analog 
computer does not subtract 20 inches from 45 inches to obtain 25 inches but, rather, it subtracts 4 volts 
from 9 volts to obtain 5 volts. This answer the operator reads as 25 inches in accordance with his ar- 
bitrarily specified "scale factor" of 1 volt equals 5 inches. In the PACE TR-10, as is the case with 
the majority of modern analog computers, the continuous variables are d-c voltages. 

The electronic analog computer makes possible the building of an electrical model of a physical system 
in which d-c voltages will behave with time in a similar way to the variables of interest in the actual 
system. If the vertical position of the center of gravity of an automobile oscillates with time during a 
disturbance, then the voltage representing the height of the center of gravity above the road surface will 
also oscillate; if the temperature of the coolant water at the exhaust port of the condenser rises exponen- 
tially to a steady value as a power plant is put into operation, then so will the voltage representing it on 
the computer. It can be said that the actual system and the electrical model are analogous in that the vari- 
ables which demonstrate their characteristics are described by relations which are mathematically equiva- 
lent. The actual system has thus been "simulated" because of the similarity of operation of the electric- 
al model and the actual system. These capabilities of analog computers are of great value in performing 
scientific research or engineering design calculations in that they give an insight into the relationship 
between the mathematical equations and the response of the physical system. Although the analog com- 
puter operates as a simulator, it performs basically as an equation solver since it performs mathematical 
operations which result in solutions of the equations used to represent the actual system. 

Although the analog computer utilizes electronic components and electrical circuit characteristics in its 
operation, it is not essential that its user have an extensive knowledge of electrical circuits. The task of 
preparing on the computer the correct electrical model is simple, and the steps necessary for accomplish- 
ing this task once the mathematical description of the primary system is known will be detailed in this 
manual. Once the electrical model is completed, experiments can be performed cheaply, quickly and with 
great flexibility to predict the behavior of the primary physical system under many different conditions, 
and it is with this in mind that one builds the model. The intention here is to give the scientist or engi- 
neer using an analog computer for the first time a comprehensive introduction to its functions, its mathe- 
matical programming, and its operation to achieve usable experimental results. 

The analog computer is basically a set of building blocks, each able to perform specific mathematical 
operations on direct voltages and capable of being easily interconnected one to another. By constructing 
an appropriately interconnected group of building blocks, an electrical model is produced in which the 
voltages at the outputs of the blocks obey the relations given in the mathematical description of a physic- 
al problem. As our interest is frequently in the dynamic behavior of physical systems, the mathematical 
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equations are usually differential equations having time as the independent variable. By applying appro- 
priate initial conditions and forcing functions to the electrical model, its behavior is determined to be the 
same as that of the primary system. This behavior is viewed through output devices which produce traces of 
voltage values plotted one against another on an x-y plotter or against time on a strip-chart recorder. Steady 
values of voltage are viewed on a voltmeter which is an integral part of the computer and can be connect- 
ed by a selector switch or by push buttons to the output terminations of many of the computing blocks. 

In reading this manual, the reader ought to gain from the first section sufficient knowledge of the com- 
puter operation to enable him to interconnect computing units and to operate the computer without difficul- 
ty, from the second section, an understanding of simple programming procedures and computing techniques 
to allow him to construct a computer circuit diagram for any straightforward investigation, and, from the 
third section, an introductory view of more involved analog computer studies which will indicate to him the 
investigations possible with a larger complement of equipment. Needless to say, the full story cannot be 
told in a short manual. There are many lessons to be learned in the art of analog computing, and, as in 
every practical subject, experience is to be preferred over reading. However, the ideas and facts pre- 
sented here ought to set the reader well along the best path to the goal of dynamic problem investigation 
by simple effective use of the modern analog computer. 
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I. THE COMPUTER 

INTRODUCTION 

The PACE TR-10 Analog Computer is a fully transistorized, general purpose analog computer. Solid state 
circuit elements are used throughout the computer to eliminate vacuum tubes, thereby achieving a compact 
design which requires very little power. It is able to operate stably and accurately in normal office sur- 
roundings. Reliable, with simplicity in functional design, it is easy to use and can be a powerful tool for 
the individual engineer in the rapid solution of scientific and engineering problems. 

Consistent with our objective of programming the computer by interconnecting its components to solve 
mathematical equations, we will concentrate our attention on the computing components — the building 
blocks for our electrical model. These blocks are able to perform the following operations on variable d-c 
voltages — multiplication by a constant, algebraic summation, integration with respect to time, multiplica- 
tien of two variables, generation of known functions of a variable, or combinations of these operations. 
Each component has input terminations and output terminations which are readily accessible at the front 
face of the computer for interconnection by plugs and patch cords. 

Below the patching area lies the monitoring and control panel. This contains features which permit (a) 
switching the computer on and off, (b) controlling the operational mode of the computer, (c) setting the 
values of problem parameters to three-place accuracy, (d) reading out stationary values of problem vari- 
ables, and (e) periodically adjusting the balance of the d-c amplifiers to ensure their accurate 
operation. 

The front face of the computer is divided into three five-inch high rows of computing components and their 
corresponding interconnecting terminations. This will be referred to as the "patch panel". In the top row 
there are attenuators for multiplying voltages by positive constants less than unity. In the bottom row there 
are high-gain d-c amplifiers uncommitted in their form of operation and, as we shall see later, capable of 
performing many tasks. In the middle row there is an assortment of components and terminations — integra- 
tor networks for use with the d-c amplifiers, fixed and variable function generators, quarter-square multi- 
pliers, comparators, and terminations for additional control panel mounted attenuators, function switches, 
the reference voltages of ±10 volts, and ground potential. These computing components will be the build- 
ing blocks in forming any model and their characteristics and uses will be explained in sufficient detail to 
enable the reader to make use of them in his problems. 

ATTENUATORS 

Probably the simplest useful operation performed on the computer is that of multiplying a variable voltage 
by a positive constant whose value lies between zero and unity -- in electrical terminology, attenuation. A 
simple potentiometer is able to perform this task effectively to an accuracy of 0.1%. Thus, one is able to 
adjust problem equation coefficients and problem initial values by setting of the appropriate attenuator. 

Four types of potentiometer modules are available, as shown in the accompanying Figure 1. 

Attenuator 42.183. Ten-turn carbon potentiometers (5K ohms) with uncalibrated knobs. The upper poten- 
tiometer has one end grounded and thus has only two terminations immediately below the knobs and on the 
left side of the module. The lower potentiometer has both ends available and thus with the arm has three 
terminations at the lower right of the module. A ground connection is conveniently placed so that when, 
as is frequently required, one end of the attenuator needs to be grounded, a short patching cable can be 
used. The arm of each potentiometer constitutes the output of the device. It is terminated a second time 
for monitoring purposes at the upper right of the adjusting knob. 
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NEAR CONTROL 
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ATTENUATOR GROUP 
TYPE 2.128 



Figure 1. Attenuator Modules 
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Attenuator 42.187. Same as for 43.183 except that in place of the second arm-termination a more convenient 
arrangement is available for adjusting the coefficient setting of the potentiometer. A push button is used 
to -replace the patched signal voltage input by +10 volts and to connect the arm of the potentiometer to a 
Pot Set Bus. With the appropriate positioning of the Meter Mode Selector Switch this connects the output 
of the attenuator directly to the voltmeter for suitable adjustment (discussed later). No auxiliary patching 
is required to adjust this attenuator accurately as is required for 42.183. 

Attenuator 42.188. Ten-turn, wire-wound potentiometers (5K ohms) with calibrated dials. Same as for 42.187 
except for a change in the potentiometer used. 

A ttenuator 2.128. This group of four ten-turn, wire-wound potentiometers (5K ohms) with calibrated dials is 
positioned at the right of the control panel. The potentiometer terminations are at the right side of the 
middle row of the patch panel. All potentiometers are grounded, having simply one input and an output ter- 
mination. Push buttons are available for ease in adjusting these attenuators and the. intention is that 
they be used in a computer model to effect changes in the important design parameters of an investigation. 

Whenever the coefficient setting of an attenuator is to be adjusted with accuracy on the computer, an elec- 
trical check of the setting must be used. A mechanical adjustment could be inaccurate for a number of 
reasons. For a fair degree of accuracy, applying +10 volts to the input and reading the output voltage on 
the voltmeter will be satisfactory, provided the intended connections from this attenuator to other compo- 
nents are complete and the potentiometer is therefore correctly loaded. For precision in setting the co- 
efficient value, a nulling arrangement must be used in which the voltage appearing at the output of the at- 
tenuator is nulled against the output voltage from an unloaded, carefully aligned, precision ten-turn poten- 
tiometer. 

The nulling arrangement for an attenuator type 42.183 is shown in Figure 2. A patch cord must be connect- 
ed from the arm of the attenuator to the VM - Jack on the control panel. A second cord must be used to con- 
nect + 10V reference voltage to the "high end" of the attenuator. With the meter mode selector switch in 
the NULL position and the reference switch associated with the NULL POT switched to the +10 position, 
the attenuator can be adjusted until a null indication is obtained on the meter. 

For the setting of attenuators Type 42.187, 42.188 and 2.128 (push button feature), a similar nulling 
arrangement is used with the* exception that the operation of the push button will automatically accom- 
plish the equivalent of the two cable connections mentioned above. The attenuator will be supplied with 
+ 10 volts, and its arm will automatically be connected to a POT BUS within the computer. Setting the 
meter mode selector switch to the POT BUS position will connect the nulling arrangement to the POT 
BUS. With the NULL POT reference switch in the +10 position, a null will be obtained by adjusting the 
attenuator. 

When setting a potentiometer it is important that the load be the same during setting as during use in the 
problem. In practice the load on an attenuator usually consists of a resistor connected between the poten- 
tiometer arm and an amplifier summing junction which rides at ground potential. In some cases it is pos- 
sible during the setting procedure to overload an amplifier so that its summing junction is no longer at 
ground. Should this happen, the overload alarm would sound and an incorrect setting would be made. To 
prevent this, the operator should ground the appropriate amplifier summing junction with a patch cord 
during the pot setting operation. The cord should be removed before proceeding with the problem. 

In developing the electrical model of the primary physical system, a computer diagram is used and there- 
fore a set of symbols is required to represent the different units. The symbol for an attenuator is shown in 
Figure 3. Note that the number of the attenuator is placed inside the circle, and the setting or adjustable 
parameter is indicated outside the circle. 
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Figure 2. Nulling Arrangement 
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a) K = adjustable parameter (setting of attenuator 
may be changed for different problem runs). 
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setting remains unchanged for a series of 

problem runs. 



Figure 3. Attenuator Symbol 



OPERATIONAL AMPLIFIERS TYPE 6.112-1. 

The most versatile computing unit within the computer is the operational amplifier. It can be used in any 
one of a number of ways to effect different mathematical operations. Sometimes, as will be seen, the opera- 
tion requires the connection of other components to the amplifier. To appreciate the use of an operational 
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amplifier consider the circuit shown in Figure 4, where the symbol containing —A represents a high-gain 
d-c amplifier. As is customary in all analog computing circuits, the voltages ei, eg, and eo are measure^ 
with respect to a common ground potential. Applying Ohm's law to the input of the amplifier, one has 
summing currents: 

H ~ e * + e ° " e s = i g (1) 
Ri Rf 



Also, by definition of A: e 0 = — A e^ (2) 

Now ig, the amplifier's input current, is extremely small compared with the currents flowing through 
Ri and Rf given by the quotients involved in equation (1). Thus, one may approximate it as zero, and, 
combining the two equations, one has: 

Rf ei 

But A, though dependent on frequency, has a value in the range of 10 4 — 10 7 for all expected computer 
operations and, therefore, provided the ratio _/ is comparatively small (less than 100, say), one can write: 

e 0 = - ei (4) 

Ri 



FEEDBACK RESISTOR 



INPUT 
RESISTOR 




Figure 4. Operational Amplifier 
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Thus one obtains the most important result that, provided the input current is low and the amplifier gain is 
large and negative, then the output- input relation is solely dependent on the ratio of the feedback to the 
input resistors. By suitable choice of resistors or other impedances such as capacitors, many worthwhile 
relations can be obtained. A few of the more commonly required ones will be detailed in the following para- 
graphs. 

1. Sign - changing or Inversion 

If both resistors are of equal value, 10K ohms say, the output voltage will be equal but opposite in sign 
to the input voltage. 

e 0 = - ej (5) 

Blue plugs containing 10K ohm resistors are available for connecting between the output and input ter- 
minations of a high-gain d-c amplifier, and single pin plugs containing 10K ohm resistors are available for 
use as input resistors. Thus the connection of a d-c high-gain amplifier for use as an inverter is as shown 
in the accompanying Figure 5. Two 100K resistors may also be used to make an inverter. However, this 
is not generally recommended since better performance is obtained with the 10K-10K inverter. 

2. Multiplication by —10. 

If the feedback resistor has ten times the value of the input resistor, then the output voltage will be op- 




CONNECTIONS PERFORMED: 
IOK 



COMPUTER CIRCUIT 
SYMBOL: 



IOK 

i 




e 0 = -ei 




-e 



Figure 5. Operational Amplifier — Inversion 
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posite in sign and ten times the input voltage. 



e a = — 10 e{ 



(6) 



Single pin 100K ohm plugs are available for use with the amplifiers and thus a connection as shown in the 
accompanying Figure 6 effects a multiplication by —10. 

3. Multiplication by a constant — a, 1 < a < 10. 

By the combined use of an attenuator and an operational amplifier, a multiplication by a constant between 
— 1 and — 10 can be obtained. Figure 7 shows the appropriate connections. 

4. Division by a constant /3. 

Division by a constant /3 can either be treated as multiplication by 1//3 or can be obtained with the cir- 
cuit shown in Figure 8. This circuit satisfies the relation: 



|8 



H_ 
Ri 



(7) 



and thus: 



1 Rf 

e o — — — — ^i 
/3 Ri 



(7 a) 



+ e ; 




AMPLIFIER 
6.112- I 
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CONNECTIONS PERFORMED: 



COMPUTER CIRCUIT SYMBOL*. 
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Figure 6. Operational Amplifier Multiplication by — 10 
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Figure 7. Multiplication by a constant — a 
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5. Summation of a Number of Voltages. 

The inclusion of more than one input resistor to a high-gain d-c amplifier circuit, each resistor having a 
voltage applied to it, changes the original development of the input -output relationship to: 

That is to say the current flowing through the feedback resistor must be the algebraic sum of currents 
flowing through the input resistors since the amplifier input voltage and current are zero. 
Thus : 

e 0 = - e x + ^ e 2 + ^ e 3 + ] (9) 

K-! K 2 ^3 

By using equal values for all resistors, one obtains a simple algebraic summation with the usu r al inversion 
associated with every computing amplifier. If the resistors have different values, then each input voltage 
is multiplied by a factor given by the ratio of the feedback resistor to the input resistor before the sum is 
taken. If both positive and negative voltages are applied to an amplifier, then due recognition to sign is 
paid in the algebraic summation process. See Figure 9. 




IOOK 



Figure 9. Operational Amplifier — Summing 
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6. Integration with respect to time. 

If the algebraic sum of the input currents is forced to pass through a feedback capacitor rather than a re- 
sistor, the current equation is: 

c tr--tk + £ + t + - 1 <10) 

This relation depends simply on the fact that when current passes through a capacitor the voltage across 
the capacitor changes with a rate proportional to the current. Integrating equation (10) and assuming an 
initial voltage across the capacitor of V 0 gives: 

If C has a value of 10 microfarads and all input resistors have values of 100K ohms-, the expression for 
e 0 reduces to: 

e 0 = - f* [ e l + e 2 + e 3 + .... ] dt + V 0 (12) 

0 

If an input resistor has a value of 10K ohms, then the corresponding input voltage is multiplied by 10 be- 
fore summation and integration. 

It is necessary to be able to control the operation of integration and also to be able to apply an initial 
charge to the capacitor to set the initial value (V o) of the output voltage eo. Thus, to form an effective in- 
tegrator more than a simple capacitive feedback is required. To effect integration on the TR-10, one 
connects an integrator network across a high gain d-c amplifier and the necessary input resistors are 
applied to that integrator network. Do not apply input resistors to the standard input terminations of the 
d-c amplifier for this would by-pass the most important "hold relay" in the integrator network (see later). 
The output voltage may be taken from either the network or the amplifier terminations. (See Figure 10). 

With these six examples there has been demonstrated the versatility of the operational amplifier. To sum- 
marize, one might list the following operations as those performed by the unit when suitably connected 
with resistors and/or capacitors. 

a) Inversion 

b) Multiplication by a constant greater 
than unity 

c) Algebraic summation 

d) Integration with respect to time of 
an algebraic summation 

The amplifier has uses other than those of directly effecting mathematical operations. By suitable analysis 
combinations of input and feedback elements can be determined for developing with one amplifier many de- 
sired transfer relations. Frequently the amplifier is used without a direct feedback connection through a 
passive element but with a feedback path through other computing components. With this arrangement the 
amplifier brings about a null of the input currents by changing that one dependent on the feedback connec- 
tion. A simple example of this use appears later in the manual when the quotient of two variables is dis- 
cussed. 

COMPUTER SYMBOLS 

Consistent with the standard symbols used in drawing computer circuits for the larger vacuum-tube com- 
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puters, the symbols used to represent summing, integrating and high -gain amplifiers are given along with 
those of other devices in Appendix I on page 67. 

CIRCUIT INFORMATION REQUIRED FOR THE APPROPRIATE OPERATION OF THE HIGH-GAIN D-C AMPLIFIER 

For accurate computation the amplifier must remain balanced. It must produce a zero output voltage when 
the combined effects of the input voltages or the absence of all input voltages demands it. With this re- 
quirement in mind the amplifiers contain chopper stabilization circuits to minimize the effects of compo- 
nent drift. Under normal circumstances the amplifiers will remain balanced for periods of hours or days. 
However, at intervals it is necessary to check this condition, and if an amplifier is found to be unbalanced 
then an adjustment is recommended. 

Below the control panel, hidden by a snap-fit cover plate, are the balance potentiometers for the amplifiers. 
With the Meter Mode Selector Switch in the position "BAL", each amplifier can be selected on the Ampli- 
fier Selector Switch and its balance tested and trimmed. When making this adjustment, it is necessary for 
all amplifiers to have a feedback connection of some form. 

In normal use the deflection of the balance meter should be adjusted to within two or three divisions of 
zero. When first turned on, the deflections may be slightly higher but will return to their normal levels 
after 30 or 60 minutes warm-up. For unusual problems which might be sensitive to amplifier unbalance or 
integrator drift, amplifiers may be balanced at more frequent intervals so as to keep the meter deflection 
below one division. 

The operation of the computer is intended to accomodate signal voltages having values up to a maximum of 
the reference voltage + 10V. Under some conditions the amplifiers can produce output voltages as high as 
± 14V, but for accurate, reliable computation they should not be required to give more than ± 10V. An over- 
load alarm feature is included in the computer which indicates the presence of an overloaded amplifier, 
usually due to a patched circuit requiring too large an output voltage. The overloaded amplifier may be lo- 
cated by placing the meter select switch in the BAL (balance) position and rotating the AMP (Amplifier 
Selector) switch until a large meter deflection is noted. The overload alarm does no more than indicate 
that something is wrong with the patched circuitry, the problem operation or the components. It is well to 
note that occasionally an overloaded amplifier will cause no appreciable error in a problem solution and 
then there is no need to take any corrective action. However, the operator should ascertain that this is 
the case before proceeding with the problem. 

MULTIPLIER, TYPE 7.045. 

Multiplication of two variable voltages is a non-linear operation which is necessary on a general 
purpose computer. A "quarter-square" technique is used to effect this operation, use being made of the 
identity: 

XY = - [ (X + Yf - (X - Yf ] (13) 
4 

The squaring operations are performed by fixed DFG's of the type described later, and interconnected 
circuitry within the unit requires only that the connections shown in figure 11 be made to the patch panel 
terminations for this unit. All four inputs ( + X, —X, + Y, —Y) must always be patched in even though one 
or both of the inputs may not change sign during a problem run. Note that the output voltage from the re- 
quired associate amplifier is — XY/10. A change in the output to +XY/10 is accomplished easily by inter- 
changing the +X and — Y inputs, or +Y and — Y inputs. 

An appropriate symbol for the multiplier is given in Appendix I, page 67. 
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DIVISION 

When it is necessary to divide one variable voltage A by a second variable voltage B, one uses a multi- 
plier in the feedback circuit of a high-gain amplifier. Consider the computer-circuitry shown in Figure 11. 
Assume the output voltage of the high-gain amplifier to be C. Then at the grid of the amplifier the null 
relationsh ip 

CB . x 

— + A = 0 (14) 

10 

is satisfied by automatic changes in the value of C. Hence, 

C = (14a) 

B 

A most important point to note about this circuitry is that although the voltage A can have both positive 
and negative values, the B voltage must always have positive values. Should B always be negative, then 
the use of — B in its place, (i.e., the interchange of the connections to +y and — y), will produce at the 
output of the high-gain amplifier, 

C = + 1^ (14b) 
B 

One must place the following restrictions on the voltages in a quotient : 

1. The absolute value of the divisor B must always be greater than or equal to the absolute value 
of the dividend A, otherwise an overload may occur in the output amplifier. 

2. The divisor B must not change sign. It must not pass through zero, for this would imply an in- 
determinate or infinite quotient. Moreover, the circuitry requires that the high-gain amplifier be 
surrounded by negative feedback, and this can only be arranged for one or other sign of the B 
voltage. If B is a negative voltage, then +B and — B are interchanged and the output of the high- 
gain amplifier is then + rather than— IS*A_ . 

B B 

X 2 DIODE FUNCTION GENERATOR, TYPE 16.101. 

As another example of the non-linear function generating equipment available in the computer, we will con- 
sider briefly the X 2 DFG. Using resistors and solid state diodes, this unit generates two quadratic curves, 
one for positive input voltages, the other for negative input voltages. Each curve is made of seven straight- 
line segments approximating the function X 2 . When used with an operational amplifier the section termin- 
ated on the upper part of the X 2 DFG module accepts a negative voltage X and produces a positive volt- 
age X 2 /10. Should X be a positive voltage, the output will be zero. The section terminated on the lower 
part of the X 2 DFG module accepts positive voltages X and produces from the associated operational amp- 
lifier a negative voltage — X 2 /10. Should X in this case be a negative voltage, the output will be zero. 

The patching interconnections required to obtain various types of mathematical functions with the X 2 DFG 
are shown in the accompanying figures 12 and 13. The symbol for use in a computer circuit diagram is 
shown in Appendix I, page 67. 

LOG DIODE FUNCTION GENERATOR, TYPE 16.126 

A second frequently occurring non-linear function which can be generated with a component available in 
the computer is the logarithm. Employing the same kind of circuitry as the X 2 DFG, the LOG DFG ac- 
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Figure 12. X 2 DFG Connections 
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cepts a voltage X with magnitude between 0.1 and 10 and with the aid of an output amplifier produces a 
voltage proportional to Log 10 [10 | X | ]. Each module contains two separate function generator circuits; 
one terminated on the upper half of the face of the module accepts negative voltages X and produces 
Y = +5 Log 10 [10 |X| ], the other terminated on the lower half of the face of the module accepts positive 
voltages X and produces Y = —5 Log 10 [lO |X| ]. For either section the absolute value of X should normally 
be restricted to lie in the range 0.1 < |X| < 10. Although input values between 0 and 0.1 may be used, 
accuracy limitations with low value inputs limit the usefulness of the answer. 

The normally required connections for the unit are shown in Figure 14. 

Just as the X 2 DFG can be used to produce its inverse, the square root, so the Log DFG can be used to 
produce the Antilog, see Figure 15* By using three Log DFG's and one amplifier, multiplication of two 
variables, neither of which change sign, is achieved. Similarly, one can obtain the quotient of two 
variables neither of which change sign. The use of one more DFG section allows either the multiplication 
or the division of a quotient by a third variable. 

VARIABLE DIODE FUNCTION GENERATOR 16.165 

In many investigations the dependence of one variable quantity Y on another quantity X is known only in 
the form of an experimentally - determined curve. For example, the drag coefficient of an airfoil is 
determined from wind tunnel experiments to be dependent on the air speed, but the dependence is known 
only as a table of corresponding values or as a graph of drag coefficient versus Mach number. The 
specific heat or thermal conductivity of a plastic material might be non-linearly dependent on the tempera- 
ture of the material and this dependence might well be known, not as a neat algebraic relationship, but as 
a graph of corresponding values. These coefficients and their dependence on the appropriate variable, 
be it velocity, temperature, or any other quantity, must be included in any computer study of the physical 
system. Their somewhat irregular dependence makes it desirable to have available in the computer a 
device which will accept a voltage representing, say, Mach number and give out a voltage representing 
the drag coefficient. For a different problem this same device will accept a voltage representing tem- 
perature and give a voltage representing thermal conductivity. The variable diode function generator is 
such a device, accepting a voltage X and when suitably adjusted, giving a voltage Y which is a required 
function of X; Y = f(X). 

Using the same techniques as the X 2 and Log DFG's, the Variable DFG permits any function of the input 
voltage X to be represented by a number of straight line segments. Solid state diodes are biased to 
change their state of conduction, each at a different value of X as X changes from —10 volts to +10 volts. 
The breakpoints (i.e. the points at which the diodes begin conducting) are at selected values of X and are 
not adjustable by the user. Where each diode conducts, it contributes a current to the output of the DFG 
unit which is proportional to the input voltage X. The constant of proportionality for each current is ad- 
justable by a "SLOPE" potentiometer associated with each diode circuit. The sum of the currents con- 
tributed by the diodes is applied to the summing junction of an operational amplifier and thereby forced to 
pass through a resistor placed around the amplifier. Thus the output voltage of the amplifier changes with 
X along a sequence of straight line segments, the slope of any segment depending on the number of diodes 
that are conducting and the settings of the slope potentiometers associated with those diodes. 

The variable DFG 16.165 has ten diodes with break points set at -9v, — 7v, -5v, -3v, -lv, +lv, +3v, +5v, 
+7v, and +9v. Each has a slope potentiometer. The value of the output voltage Y = f(X) at X = 0 can be set 
at any value between — lOv and +10v by the PARALLAX potentiometer, and the slope of the curve at X = 0 
can be set by the CENTRAL SLOPE potentiometer to have any value between plus and minus 2 volts per 

* Amplifier output signals are obtained by summing Amplifier input currents. The current output of the 16.126 Log 
DFG is 10 " 3 log 10 (10X) amperes. 
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Figure 17. Nonlinear Function Easily Generated by a Variable DFG 



volt. The slope potentiometer of each diode allows the slopes of adjacent straight line segments to differ 
by any value up to 1 volt per volt. 

For the generation of any arbitrary function, two uncommitted operational amplifiers are required to be 
connected to the variable DFG as shown in Figure 16. This arrangement would permit the function shown 
in Figure 17 to be generated quite easily. 

The adjustment of the DFG to produce a desired function is facilitated by the use of an Attenuator unit 
type 42.243, which by connection to an inverting amplifier produces a voltage which can be accurately 
stepped from —10 volts to +10 volts in 1 volt increments. This permits the input voltage to be set con- 
veniently to any of the breakpoint values, and then the slope potentiometers can be adjusted until the 
amplifier connected to the output of the VDFG reaches the desired voltage f (X). If no Attenuator unit is 
available the required breakpoint voltages have to be obtained from one of the common attenuators con- 
nected to an inverting amplifier. 
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The following steps should be followed in preparing a variable DFG to give a desired function (see figures 
18 and 19). 

1. Prepare a table of the appropriately scaled voltages desired at the output of the DFG when the 
input voltage has each of the following values: —10, —9, —7, —5, —3, —1, 0, +1, +3, +5, +7, +9 and +10 
volts. Note that the values at —1, 0 and +1 must be colinear. 

2. Plug the VDFG into a service shelf and connect the required amplifiers appropriately to the ter- 
minations of the unit. The VDFG potentiometers are now available for screw driver adjustment. Turn 
all potentiometers fully counter-clockwise. 

3. Ground the VDFG input or set the input to zero and, by adjusting the PARALLAX potentiometer, 
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Figure 18. Patch Panel Arrangement for Adjusting a Variable DFG 
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Figure 19. Circuit for adjusting a Variable DFG 



obtain the appropriate output voltage f (0) from amplifier 2. Use the procedure detailed later on page 32 
to check this voltage. 

4. Apply —1 volt to the VDFG input and- adjust the potentiometer labeled "—1" until the output of 
amplifier 2 reaches the value f (— 1). 

5. Apply —3 volts to the VDFG input and adjust the potentiometer labeled "—3" until the output of 
amplifier 2 reaches a value f (— 3). 

6. Continue to change the input voltages and by adjusting potentiometers obtain the appropriate 
values at the output of amplifier 2 as listed in the following table: 
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7. Quickly check through all points to ensure that the complete function is well adjusted. 

8. Remove the service shelf and insert the VDFG into its operating position. 

It is well to note that for certain functions which are monotonic and for which the slope is always in- 
creasing with the absolute value of X { e.g. X 3 , e x , Tan X }, the VDFG requires only one amplifier. The 
appropriate connections for this situation are shown in Figure 20. 

VARIABLE DIODE FUNCTION GENERATORS 16.154 AND 16.156 

The variable DFG 16.154 accepts only negative input voltages (— 10 < X <0). The variable DFG 16.156 
accepts only positive input voltages (0 < X < +10). Each has nine diodes with breakpoints set at integral 
values of voltage 1 through 9. Its operation is quite similar to that of the type 16.165 and the required 
connections and adjustment procedures correspond closely. In the adjustment procedure the PARALLAX 
and ZERO SLOPE settings are made first, before proceeding to the diode SLOPE settings which are made 
in order. When obtaining monotonic functions the "one amplifier" connection shown in Figure 20 may 
also be used with VDFG's 16.154 and* 16.156. 

In some cases it may prove desirable to generate a function over the range — 10 < X < +10 with more ac- 
curacy than can be obtained with VDFG 16.165. By connecting a 16.154 and 16.156 in parallel the 
necessary function can be generated with 19 segments rather than with the 11 segments available with 
16.165. 

MISCELLANEOUS DEVICES AND PATCH PANEL TERMINATIONS. 

In addition to the standard computing components, a number of devices are found to be useful in program- 
ming a problem investigation on the computer. They can be used by the operator as he desires and are re- 
viewed here to complete the list of items terminated at the patch panel. 

1. Dual Function Switch, Type 2.127. 

Switches mounted on the control panel are terminated at points in the middle row of the patch panel. Al- 
most self-explanatory in their use, these function switches are single-pole, double-throw with an OFF 
(disconnected) center position. 

2. Signal Voltage Comparator, Type 6.143. 

An automatic switching device, the comparator operates two sets of contacts, permitting computed prob- 
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Figure 20. Variable DFG Connections for Monotonic Functions 



lem voltages to determine connections or conditions applying in a patched circuit. As its name implies, 
the comparator accepts two input voltages, compares their sum to zero (approximately) and positions two 
switches up or down, depending on whether the sum is greater than or less than zero. The details of the 
unit's terminations are given in Figure 21. 

An appropriate symbol to represent the comparator in circuit diagrams is given in Appendix I, page 67. 
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4. Adjust this potentiometer until relay switching occurs. 
Bias (IN 2 ) input is now set for proper switching level. 

5. Remove coefficient potentiometer from IN x and patch 
variable X into IN X . 

6. The relay contacts will be in the positive position when 
X > K, they will be in the negative position when 
X<K. 



Figure 21. Signal Voltage Comparator 
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3. Overload Alarm, Type 13.012. 

This unit provides an audible warning signal of about 400 cps when an overload occurs in any of the 
operational amplifiers, i.e. when the summing junction error exceeds a safe level for any cause. When the 
alarm operates, the overloaded amplifier can be located quickly by setting the meter selector switch to the 
BAL (balance) position and rotating the AMP (Amplifier Selector) switch until a large meter deflection is 
observed. 

4. Reference Voltage Supply . 

The computer uses plus and minus 10 volts as sources for all computer signal voltages. Balanced about 
ground potential, these sources are available at the patch panel and are also connected internally to many 
points throughout the computer. 

CONTROL PANEL 20.344. 

In order to allow simple control and monitoring of the computing components, this panel contains the fol- 
lowing components : 

1. Primary power switch. Positions : ON, OFF. 

2. Neon light. This light works in conjunction with the primary power switch and indicates when 
the computer is switched on. 

3. Mode Control Switch. Positions : RESET, HOLD, OPERATE. 

In RESET any circuit patched on the computer is functioning except that the outputs of all integrators are 
held at their required initial conditions. The programmed problem is thus set to those conditions corres- 
ponding to time zero. 

By switching to OPERATE, the integrators are simultaneously set free and with the voltages applied to 
their input resistors causing changes in the output voltages, a time-varying behavior is produced. This will 
be the voltage solution of the programmed problem. 

Switching to HOLD permits the solution to be held or "frozen" at any time that is convenient. After making 
whatever observations that dire desired, the solution can be resumed by switching back to OPERATE or re- 
turned to its initial conditions by switching to RESET. 

4. Voltmeter. Used as a voltmeter, null meter or balance indicator. 

5. Meter Range Switch. This switch permits the selection of the following meter ranges when the 
meter is used as a voltmeter (i.e., in the AMP or V.M. positions of the Meter Mode Selector Switch) : 

± 30 volts 
+ 10 volts 
± 3 volts 
+ 1 volt 
± 0.3 volts 
±0.1 volts 

6. Meter Mode Selector Switch. Positions : POT BUS, NULL, V.M., AMP, BAL. 

POT BUS — Connects the meter as a Null meter to the computer Pot Bus for setting attenuators 42.187, 
42.188 and 2.128 (push button feature). 

NULL — Connects the meter as a Null meter to the VM jack connection on the control panel. This po- 
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sition is used for setting attenuators 42.183, which do not have the push button feature, and for measuring 
accurately other output voltages in the computer. 

V.M. — Connects the meter as a voltmeter to the VM jack. Voltages patched to the VM jack may thus be 
read on the meter. Different meter ranges may be obtained with the meter range switch. 

AMP - Connects the meter as a voltmeter to the AMPLIFIER SELECTOR SWITCH so that the output vol- 
tage of each amplifier can be read on the meter. 

BAL — Connects the stabilizer outputs of the operational amplifiers, as selected by the AMPLIFIER SE- 
LECTOR SWITCH to the meter, which is now connected as a high sensitivity voltmeter. Thus, in this po- 
sition, the stabilizer outputs of all amplifiers may be read on the meter and balanced to zero, through the 
use of the amplifier balance adjustments. 

7. NULL POT and Reference Switch. Used to apply a precision voltage to the Null meter for set- 
ting attenuators or reading output voltages accurately. The reference switch controls the polarity of the 
voltage applied to the NULL POT (±10 volts). The switch should be in the positive position for nulling 
positive voltages or setting attenuators and the negative position for nulling negative voltages. The NULL 
POT and reference switch operate in conjunction with the POT BUS and NULL positions of the meter se- 
lector switch. 

8. Amplifier Selector Switch (AMP). The first 20 positions of the switch permit the selection of the 
outputs of the operational amplifiers. The voltages at these outputs may be read on the meter by placing 
the METER MODE SELECTOR switch in the AMP position. A second deck on the amplifier selector 
switch allows the selection of the amplifier stabilizer outputs which may be balanced to zero with the 
METER MODE SELECTOR switch in the BAL position. 

Switch position 21 selects the positive reference voltage (+10V). Switch position 22 selects the negative 
reference voltage ( — 10V). 

9. VM Jack. This terminal is connected to the meter when the METER MODE SELECTOR switch 
is in the NULL or VM position. 

10. AMP OUT Jack. Amplifier outputs, as selected by the AMPLIFIER SELECTOR SWITCH, are 
available at the jack at all times for metering or recording purposes. 

As pointed out in Paragraph 6 above, any steady voltage existing in the computer circuit may be read to 
three place accuracy by using the NULL POT in conjunction with the NULL position of the Meter Mode 
Selector Switch. 

EXAMPLE: To read the output of an amplifier to three place accuracy : 

a. Connect AMP OUT Jack by means of a patch cord to VM Jack. 

b. Select amplifier output with the AMPLIFIER SELECTOR switch. 

c. Switch METER MODE SELECTOR switch to the NULL position. 

d. Observe deflection of meter. If needle deflects to the right (positive), switch null pot reference 
switch to +10. For a negative deflection switch null pot reference switch to —10. 

e. Adjust NULL POT until a null indication is obtained on the meter. 

/. To obtain the value of the amplifier's output voltage, multiply the NULL POT dial reading by 
the reference voltage used in step d. 
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REPETITIVE OPERATION 

Repetitive Operation is an available feature which is most useful when solving certain kinds of problems. 
The problem solution time is reduced by a ratio of 100 : 1 and the computer automatically alternates be- 
tween the RESET and OPERATE modes, causing the time behavior of interest to be produced many times 
each second. A solution which is programmed under normal operation to be obtained in 5 seconds, is 
achieved in 50 milliseconds, permitting it to be repeated approximately twenty times each second. The 
high-speed solution can be viewed on an oscilloscope and adjustments made to the potentiometer set- 
tings, their effect being seen immediately. Quickly a desired behavior can be sought, and then by return- 
ing the computer to normal operation it can be permanently recorded at slower speed on a plotting table. 
The time taken for one solution to be displayed is so short compared with the standard human recognition 




a) Control Unit Type 20.532 b) Integrator N etwork Arrangement for Repetitive Operation 

Figure 22. 
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time that changes in behavior from solution to solution appear to be continuous; the behavior appears to 
sweep towards that desired. 

To achieve repetitive operation three units are required. 

1. A control unit type 20.532, which fits into the computer's control panel. It contains a multi- 
position switch, which switches the computer over to repetitive operation and also determines the total 
time of solution — 10, 20, 50 and 100 millisecs. A vernier adjustment permits the time of solution to be 
further adjusted between these values; with the multi-position switch turned to 20, say, moving the vernier 
knob from its fully counter-clockwise position changes the time of solution gradually from 20 to more than 
50 milliseconds. The control unit has a termination at which a sweep voltage for any standard oscillo- 
scope is available. The use of this voltage eliminates difficulties in synchronization that would other- 
wise need to be overcome. (See Figure 22a.) 

2. A timing unit type 36.082 which fits in the back of the console next to the reference power supply, 
and supplies the alternating, variable duration, relay control signals to the integrator networks. 

3. A set of integrator networks type 12.425 which in place of the integrator networks type 12.263, 
contain both 10 microfarad and 0.1 microfarad capacitors permitting the 100:1 change in speed of solution. 
These units have the same terminations as the type 12.263, illustrated earlier. With the multi-position 
switch of the control unit switched off the 10 microfarad capacitors are in the circuit. With the multi- 
position switch in any other position, the 0.1 microfarad capacitor is in the circuit. (See Figure 22b.) 

In place of the standard "HOLD" and "RESET" relays shown in Figure 10, a high-speed electro-mechani- 
cal relay, controlled from the timing unit, switches the summing junction between the mid-point of the 10 K 
ohms initial condition resistors (RESET) and the input resistors (OPERATE). The integrators are RESET 
for 10 milliseconds between each solution. The time for which the integrators OPERATE is controllable 
between 10 milliseconds and approximately 250 milliseconds. During the setup of a problem when at- 
tenuators are being adjusted the Rep Op control switch should be turned to OFF so that internal relays 
will ground the summing junctions. If this is not done, loading errors will cause the attenuators to be set 
incorrectly. 
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II. BASIC PROGRAMMING 

The earlier sections of this manual emphasized the use of the basic computing elements of the TR-10. 

An integrated review of the steps necessary to program a problem for analog computation will now be dis- 
cussed. 

The general procedure is : 

1. Obtain mathematical equations for the system to be simulated. 

2. Bring the equations into a convenient form for solution on an analog computer. 

3. Develop a simplified block diagram to indicate the basic signal flow in the computer program. 

4. Perform amplitude and time scaling on the original equations. 

5. Develop final computer diagram. 

6. Perform static check calculations. 

7. Patch the problem, monitor static check voltages, and run the problem. 
SYSTEM EQUATIONS 

Starting from the physical system one can derive a mathematical description using the commonly-known 
physical laws, experimental results and any reasonable simplifications that might appear desirable. The 
description will include a set of equations, coefficients values, and possibly some experimental curves 
relating variables in the study. 

Consider now the available analog computing equipment and decide what equation modifications, if any, 
are necessary (for instance, linearizing wherever possible to reduce the complexity of the study, using 
finite differences to reduce the number of continuous independent variables in a partial differential equa- 
tion). Make any changes of variables or formal modifications that facilitate programming. Remember that 
the equations most conveniently solved on the computer are ordinary differential equations. 

REARRANGING THE EQUATIONS 

The system equations (usually differential equations) must now be brought to a form in which an analog 
computer of the general type can solve them. The original equations are rearranged to have only the term 
with the highest derivative on the left-hand side. 



dt n X dt n ' 1 dt n ' 2 



(15) 



Successive integration of the highest derivative will yield all the lower derivatives on the right-hand side 

of the equation. For solution on the computer it is assumed that ^ 7 exists as' the input to an integrator; 

dt n 

by producing the lower derivatives by integration one can sum appropriately all terms on the right-hand 
side of the equation to form d n y . This sum of. terms must therefore be fed into the integrator where d n y 
was assumed to exist. u ' " 

Example : 

an a mi ati r\rt • u . { 



Given equation : d 2 y ^ dy v i f \ Qg) 



Rearranged equation: d 2 y _ a 2 dy a 3 + a * F (t) (17) 

dt 2 dt a± a x 
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Assume ^ exists and integrate twice to form ^Xand y. Generate Fit) and form the R.H.S. (right-hand 

side of the equation (17) by summation. Use this sum for ^L. . 

dt 1 

BLOCK DIAGRAM 

The block diagram is a preliminary computer diagram ignoring scales, signs, gains, but showing computing 
elements and their interconnections. It is frequently useful for obtaining an estimate of" the amount of com- 
puting equipment needed, and for producing at a later time a neat computer circuit diagram. 



F(t) 



d 2 Y 



I 



dY 
dT 



/ 



Figure 23. Block Diagram for Equation 17 



SCALING 

Depending on the dimensional units used, variables can have values or rates of change that range from 
very small to very large numbers. When preparing a problem for solution on the computer, it is necessary 
to arrange that voltages representing these variables will never have values larger than the allowable 
value of ± 10 volts, nor will they change rapidly enough to exceed the frequency limitations of the com- 
puting and recording equipment. 

The maximum values which quantities can assume in a particular physical system can usually be esti- 
mated from i 

a. A knowledge of the physical system under consideration. 

b. An investigation of the system equations. Steady state relations and characteristic frequencies 
may be employed to evaluate the ranges of the dependent variables. 

AMPLITUDE SCALING 

For maximum accuracy, the signal voltages representing the dependent variables and their derivatives 
should be as large as possible without causing the equipment to overload. The equipment is designed for 
signal voltages within the range of + 10 volts. The maximum values of the physical variables must there- 
fore be scaled so that: 

(maximum value of physical- variable) x a = 10 volts (18) 

whereat ™ volts Q9) 

max. value of physical variable physical unit 

is called the amplitude scale factor. 
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The following examples demonstrate how to determine amplitude scale factors for a problem : 

Physical Variable Scale Factor a Scaled Variable 

Linear velocity x ^ 100 ft/sec JL volts r^i 

10 ft /sec L 10 J 



Angular acceleration 0 = 1 JL rad/sec 2 



$ i Km ft 1 volts 

t /sec 

volts [5 e ] 



Linear displacement x = 1500 ft J_ volts r x i 

200 ft /sec 200 



2 rad/ 



sec" 



Temperature T = 800° F — [— ] 

100 °F 100 



To identify scaled variables which wil| appear as voltages in the computer circuit, square brackets shall 
be used. 

TIME SCALING 

The ideal duration of a single run on the computer is somewhere from 15 to 60 seconds. Too long a solu- 
tion time is inefficient and may introduce small errors due to integrator drift. The minimum time for a com- 
puter run is normally determined by the speed of response of the computing and recording units. If mechanic- 
ally driven recording equipment is used, the frequencies occurring in the problem setup should be of the 
order of 1 to 2 cps. To satisfy these recommendations a change in time scale can be introduced by apply- 
ing a time scale factor /3, so that 

r = fit (20) 

t = machine time in seconds 
t = real time in seconds 

Furthermore, dx dx dr dx . . 

= = p v^l j 

dt dr dt dr 

Let us again consider an equation of the form : 

P- = - - K 2 y + K 3 F(t) (22) 

dt 2 dt 



If we set = z we obtain 
dt 



P = - K t * - K 2 y + K 3 F(t) (23) 
dt 



Furthermore let it be assumed that the amplitudes of the variables were scaled as follows : 



[ «i — ] = [ a 1 — ] ; [ a 2 -2-] = [ a 2 z ] , [ a 3 y ] , [ a 4 F ( t) ] 
dt 2 dt dt 



= [«, -1; la. 4- 

dt 2 dt dt 

Taking account of the amplitude scale factors and time scaling the equations to be programmed on the com- 
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puter become 



BU^-] [a 2 z] [a 3 y] U 4 F(^-)] 
= — K x — K 7 + K 3 



a 3 



a-4 



(24) 



or 



r „ <iz i 1 oti^i 1 a 1 K 2 _ 1 a-,K 3 r / r s 

L Gl — J C a 2 z ] - - -i-?- [ a 3 y ] + _J_i [ a 4 F( — ) ] 

dr P a 2 P a 3 Z 3 a 4 £ 



(25a) 



and 



[ a 2 — ] = — . [ Oj2r ] 



(25b) 



dr 13 
COMPUTER DIAGRAM 

The resulting computer circuit, appropriately scaled, is shown in Figure 24. 

Note in the circuit shown in Figure 24 that every input to an integrator is divided by the time scale factor 
(3 . It is possible, in fact, to establish the following general rule : 

For a change in time scale every input to an integrator must be multiplied by i. 

Jl < 1 for slowing down a problem. 

.L > 1 for speeding up a problem. 

P 

STATIC CHECK CALCULATIONS 

The static check is a procedure ensuring that the problem set on the machine corresponds to the original 
equations. A set of arbitrary check values are chosen for the variables and their (n-l) derivatives and 
the value of the n 1 * 1 derivative (left side of the rearranged equation) is calculated. When the problem is set 
up on the computer, initial conditions corresponding to these check values are imposed, the voltage for the 



0C4 



-[«c 4 r«%3— O 



_cy q_ 



OC.K 



ri 



.JSjJUQc.y,,,] 



-hH-i[«'-H-] 




O 



GC, K 



1 "2 



<*3 



(OC3Y (T)] 



Figure 24. Computer Circuit Representing Equation 25 



38 



TR - 10 OPERATOR'S MANUAL 



n { h derivative is measured and compared with the previously calculated scaled value. Various other vol- 
tages appearing at outputs of amplifiers should also be calculated and measured. 

EXAMPLE OF AMPLITUDE AND TIME SCALING 

a. Scaling Initial Conditions. 

If an integrator delivers, say — [5*] and the initial condition is x (t - 0) = X 0 = 1.2, this is realized by 
introducing a constant voltage at the IC - terminal of the integrator, taking account of the sign reversal 
between IC input and integrator output. 

The constant voltage is obtained from a potentiometer fed by ± 10 volts reference. The output of the po- 
tentiometer should be (5) {X 0 ) = (5) (1.2) = 6 volts, and the potentiometer setting, considering the 10 
volt input, is therefore 

5x 0 



10 



= 0.600 



(26) 



Since x 0 is positive and the output of the integrator is defined as — [Sx] , we must take the + reference 
voltage. 




X (t= 0) =X Q = 1.2 



[5X] = - [5X Q ]= -6 VOLTS AT t = 0 



Figure 25. Scaling Initial Conditions 



b. Scaling Variable Inputs . 

If it is required to generate a function 

J = g(t) (27) 
as an input into the computer system, we scale the function as follows : 

First we must find the maximum value of y and amplitude scale accordingly, say [20y]. The time scale 
relation is r = (3 1 . Therefore: 

[20y] = 20g (^) = gi (r) (28) 

The computer circuit or external equipment producing g(t) must be arranged to deliver after r seconds 
running time an output of g 1 (r) volts. 



c. Explicit Time Dependence 

If the original system is of the form 



dXj 
dt 



= fi (Xi ... X n , t) 



(29) 



39 



TR-10 OPERATOR'S MANUAL 



GIVEN DATA 
-O.I8ft. < X 



MAX 



< 0.18ft 



K s 420 LBS/FT 

m = 0.08 SLUGS 
LBS 



C = 2.0 



FT/SEC 



Figure 27. Spring-mass-damping System 



d 2 x 
dt 2 



'm 



dx 
dt 



m 



+ X 



Figure 28. Block Diagram for Equation 33 



natural frequency of the system, 

- ■ ■ il - 73 (34) 

Assuming a solution of re = sin co t (35a) 

we obtain: dx , . . 

— = a> cos ot (35b) 
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d 2 x 
dt 2 



= — co 2 sin cot 



(35c) 



Thus, for a maximum amplitude x max = 0.18 ft we can estimate for a damped system the maximum values 
for the velocity 



and the acceleration 



^ < 0.18 co = (0.18) (73) « 13 ft/sec 
at 



— < 0.18 co 2 = (0.18) (73) 2 « 950 ft/sec 2 
dt 2 



To represent these variables by voltages which may not exceed ± 10 volts we have to apply amplitude 
scale factors so that the scaled variables become: 

[50,]; [£]; [^] 

TIME SCALING 

From the original equation we find that the maximum system frequency can be approximately 12 cps. In 
view of the recording equipment used a time scale of jS = 10 seems a good choice. Every input to an in- 
tegrator must therefore be multiplied by = -~ which will result in a computer solution which is ten 

times slower than the real physical process. 

We can now write the scaled equation by simply introducing the amplitude-scaled variables into the origin- 
al equation. At the same time we shall make proper adjustments to keep the equation true. We may neg- 
lect the time scale factors /3 in the process of writing scaled equations, as long as we remember to add 
it in the computer circuit. 



Scaled equation 



[-*-] = -L- | -£(2) [^] Ll^ll (36a) 

100 . 100 I m 2 m (50) J 



[_*_] = _ 1<L [*] K— [S0x] (36b) 

100 100m 2 5000m 



FINAL COMPUTER DIAGRAM: 

From the scaled equation and with the help of the block diagram we can now develop the final computer 
circuit shown in Figure 29. 

Note that the settings of all potentiometers in the above circuit are less than unity. This is accomplished 
by dividing the pot settings dictated by the scaled equation by an appropriate constant. To maintain the 
equation these constants were compensated for by using an amplifier gain of the same magnitude. 

STATIC CHECK CALCULATIONS 
The original system equation is: 



d 2 x 
dt 2 



m dt m 



(37) 
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+ 10 



For static 
check only 



- 10 

50 Xp 
W TO" 




5000m (10) 

Figure 29. Computer Circuit for Spring-mass- damping System 

Assume for the purposes of a static check that: — 8; x = 0.1 

at 



d 2 x 
dt 2 



2 ( 8 ) _ (o.l) = 725 



0.08 



0.08 



(38) 



If the problem was correctly programmed and patched, the following voltages should exist in the computer 
circuit: 

x - 725 

+ [ — ] = = - 7.25 volts output of amplifier 1 

100 100 

x 8 

— [y] = — -g = - 4.00 volts initial condition on amplifier 2 

+ [ 50%] = (50) (0.1) = +5.00 volts initial condition on amplifier 3 

+ [ = +4.00 volts output of amplifier 4 

2C x 

[— ] = +2.00 volts output of potentiometer 21 

100m 2 



K 



5 x 10 4 i 



[50*] = +0.525 volts output of potentiometer 22. 



The Voltmeter may be used to monitor these voltages. After the static check has been completed the com- 
puter is ready for operation to investigate the behavior of the mass-spring system under various initial 
conditions or changes in the values of m, C and K. 

SOLUTION OF SIMULTANEOUS DIFFERENTIAL EQUATIONS 

The programming of a set of simultaneous differential equations presents no more basic difficulties than 



43 



TR-10 OPERATOR'S MANUAL 



the solution of a single equation. The following problem serves as an illustration. 
AN AUTOMOBILE SUSPENSION SYSTEM 

Assume the task is to investigate the response of an automobile suspension system for selected disturb- 
ances. The system response to these disturbances for various values of the system design parameters 
shall be obtained. 



Physical Constants 

M l = 25 slugs 

M 2 = 2 slugs 

K, = 1000 lb/ft 

K 2 = 4500 lb/ft 

D = Variable 



Ranges of Variables 



^ 6 in. 
x 2 = 12 in. 
x 3 = 6 in. 

(See Figure 30.) 



x 3 =f(t) 




777777 



M 2 

K x 

K 2 

D 

x l 

x 2 

x 3 



One quarter of mass of automobile 
Mass of one wheel and half of the axle 
Constant of main auto spring 
Spring constant of tire (assumed linear) 
Shock absorber damping constant 
Displacement of auto body 
Displacement of wheel 
Roadway profile displacement 



Figure 30. Simplified Automobile Suspension System 



SYSTEM EQUATIONS 

The differential equations of motion of the system are derived by equating forces acting upon the masses. 
Thus for M t : 

Mi d^_ + D ( dx>_ _ dx^ + (%i _ x j = 0 (39) 
dt 2 dt dt 



and for M 2 : 



M 2 #*2 + D { ^2_ _ dx^ + Ri {Xi _ Xi) + Ki {x2 _ Xi) = o 
dt 2 dt dt 



(40) 



Estimate of maximum values: 



The undamped natural frequency of Mj, is 
OJ ni - v — - = 



V looo 

25 



= 6.32 rad/sec 



(41) 
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For the mass M 2 we obtain : 



con* = V Kl * K > = V 5500 = 52.4 rad/sec (42) 
Reasonable estimates for the maximum values of the variables and their derivatives are: 



Ximax 


= 0^5 ft 


Ximax 


- \ ft x 3T nax = 1 ft 


x \max 


~ ( x imax ) 


Oinx = 


(0.5) (6.32) « 3.2 ft/sec 


• 

x imax 


~ ( x imax ) 


0t) w2 = 


(DC 52.4) « 50 ft/ sec 


x imax 


= (*imox) 


<»ni K 


20 ft/ sec 7 


x imax 


= (^jmox) 




2590 ft/ sec* 



We can therefore scale the variables so that the following voltages will appear in the computer circuit: 

[20*J [10 *J [10*,] 

[2*J [^-] 
5 

2 500 

Solving the original equations for the highest derivatives and introducing the scaled variables yields the 
scaled equations: 

r -. n 5D r X 2 -. D r . , &i r n ^1 r 

[ Y ] " 2M7 [ ^ ] " SfT [ & ' 1 + 20MT [ 10 *- 1 - 4ST [2(k ' ] 

r i 2 , D r . , D r Xo Ki r , K, + K 2 r , K 2 r -| 

[ = ____ [ & ] _ [ — ] + — [ 20x 1 ] - — [ 10* 2 ] + — [ 10*3 1 

500 1000M 2 100M 2 5 10«M 2 5000 M 2 5000M 2 ^ 

The coefficients in these equations can each be separated into two parts, — an attenuator setting and an 
amplifier gain. Assuming a range of values of D up to 100 we have on rewriting the equations: 

[^] = 10 (— )[^] - (— )[2*J + 10(-^-)[10* 2 ] - (-A-)[20xJ 
2 4M X 5 4¥i 200^ 40 M x 

rA i _ ( _ n <£ ] _ ( + ( _^_)[20* 1 ] - (^-^)[10*J + ( -^-)[10* 3 ] 

1 500 J " ^ 1000M 2 H 1 ^ 100M 2 H 5 ^ 10<M 2 5000 M 2 5000M 2 

(44) 

From the scaled equations we can now develop the computer circuit shown in Figure 31. An attenuator 
placed immediately before an integrator must produce with the gain shown on that integrator the 

necessary change in scale factor across the integrator, e.g. a setting of % followed by a gain 10 causes 

x ' 1 ' 

— [ -g-] to be integrated to + [10# 2 ]. The factor — permits the inclusion of a time scale change. In this 

problem /3 = 10 appears appropriate, slowing down the solution by a factor of 10 to 1. 

The disturbances x 3 = f(t) may be fed to the computer from a recording of road noise, or the voltage [10* 3 ] 
may be a suitable artificial disturbance produced on the computer -- a step, ramp or sinusoidal time 
function. 
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The voltages representing x x and x 2 can be recorded on a strip -chart recorder or on a VARIPLOTTER. 
Typical curves of time behavior for an input x 3 chosen as a 35 millisecond pulse are shown in Figure 
32. The effect of changing the value of D is obvious. 

A NON-LINEAR PROBLEM 

The problem is to compute a curve of radial velocity of a particle vs radial distance from the center of a 
disc used in a spray drier. The results of the computation shall be used to select the design parameters 
of a drier for a particular duty. 

The differential equation which describes the velocity of the particle is : 

V — + AV 3 - Br = 0 (45) 
dr 

V = Radial velocity [in/sec] 

r = Radial distance from center of disc [inches] 

A, B, = Constants which in the case of laminar flow depend upon the angular 
velocity of the disc, volumetric feed rate, vane height and the density 
and viscosity of the liquid. For turbulent flow the same equation holds, 
with different values for constant A. 

The velocity will be recorded as a function of the radius r on a VARIPLOTTER. 

Considering the original equation we must replace the independent variable r by time t, since the com- 
puter integrates with respect to time only. Let us therefore introduce 

r = kt (46) 

For k = 1 we obtain : 

£ . *i _ ^ <47) 

dt V 

The constants A and B will be varied for several runs on the computer. Their mean values are : 

A = 1.66 x 10 " 4 
B = 3.6 x 10 7 

The velocity V will have a range of : 

V = 0 6000 in/sec 

for a range of radius r = 0 .... 5 inches. We can therefore scale the variables as : 

[2t]; [10-* V] 

With the scaled variables we obtain the scaled equation as follows: 

Introducing amplitude scale factors and maintaining the terms of the equation consistent: 

[10-3^] = 121 B r i2t \ - — [10-K]» (48) 

Including the factor of ten which occurs in the output when doing multiplication and division (see page 15), 

[10-3 dV] _ IP" 7 B [10] [2t] A_ [10- 3 F] 2 (49) 

dt 2 [10-»K] 10- 4 [10] 
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Dividing the coefficients where necessary into two parts, one less than unity, the other a power of 10, 
so that they can be obtained using attenuators and amplifier gains, 

nn-s rfF l in t B \ [ 1Q ]I>] nmnn [10-' VV (cm 

[1 ° H ] = 10 ( 2Tib". ) TT^TT " 10 (1000A) lIor (50) 

From this scaled equation we see that we have to generate a voltage equal to 2 t. This is achieved by 
feeding an integrator with a voltage of constant value — 2 volts. 

The computer circuit can now be developed from the scaled equation and is shown in Figure 33. 

This circuit has no unusually large amplifier gains nor unusually small attenuator settings. Thus a time 
scale change appears unnecessary. Indeed the full solution time, as we know already, is five seconds, 
and this might be perfectly reasonable and not give rates of voltage change within the circuit which are 
inappropriate for the plotting equipment to be used. If a slower solution is found to be desirable, then it 
can be achieved by the accepted method of reducing the gains associated with each integration. The gain 
change is made either by resetting a potentiometer or changing the input resistor to the amplifier, which- 
ever is the more convenient. For example, to double the time of solution the settings of attenuators 3, 4, 
and 8 can all be divided by /3 = 2. 



IOOOA 




Figure 33. Computer Circuit Representing Equation 50 



49 



TR-10 OPERATOR'S MANUAL 



To convert the voltage at the output of integrator 2 into inches we must remember that t — r and thus 2 
volts represents 1 inch. The maximum output of the integrator will be 10 volts for r = 5 inches. The 
maximum output of the velocity integrator will be 6 volts for V = 6000 in/sec. We shall therefore use an 
arm scale (x) on the VARIPLOTTER of 2 volts/inch, while the pen scale (y) will be 1 volt/inch. 

The calculation of the static check voltages should be performed similarly to the procedure indicated 
previously. 

Typical curves obtained from the computer for several values of B with A held constant and vice versa 
are shown in Figure 34. 




2 3 4 

RADIUS IN INCHES 
(o) 




2 3 4 

RADIUS IN INCHES 
(b) 



Figure 34. Plot of Radial Velocity as a Function of Drier Radius 



GENERATION OF ANALYTIC FUNCTIONS 

The following is a summary of some typical computer circuits for the generation of analytic functions 
sometimes required in the investigation of physical problems. 

■Xt 



1. Generation of y = e 



— - — Ae^ f = — Ay 
dt 



2. Generation of sin co t and cos co t 

d_ 
dt 

d 



sin co t = co cos co t 



_ cos cot = — co sin co t 
dt 

sin co t = co f cos cotdt 



(See Figure 35.) 



(51) 
(52) 

(53a) 
(53b) 
(53c) 
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+ 10 




Figure 35. Generation of y - C" 2 1 



cos cot = — co f sin cot dt 

(See Figure 36.) 

3. Generation of y = e~ x where x and y are both dependent problem variables. 

dy _ dx -x 



dt 



dt 



dx 



Here it is assumed that -gj- = x is available as the input function to this computer circuit. 

(See Figure 37.) 

4. Generation of y = In (x + a) 



dt x + a dt 



(See Figure 38.) 



(53d) 



(54) 



(55) 



CURVE FOLLOWER 2.040 



The VARIPLOTTER, with the addition of a curve follower head, may be used to produce an arbitrary or 
non-analytic function of an input variable. The function is drawn on graph paper with a special con- 
ducting ink. One end of the* curve is then connected to the output of a 470 kc oscillator and the other 
end is grounded. The curve (conductor) will generate a magnetic field which on one side of the curve is 



-10 



Cos(O) 



= O.I 



+Cos tut 




Figure 36. Generation of sin co t and cos co t 
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180° out of phase with the field on the other side of the curve. (A curve follower head is mounted in 
place of the pen of the VARIPLOTTER.) The curve follower head will sense an excursion from the 
curve and will supply an error voltage to the pen servo-system in the VARIPLOTTER, which will return 
the curve follower to the center of the line of the curve. The arm of the plotter is driven by an external 
voltage x. The output voltage corresponding to y = f(x) is obtained from the wiper of the pen follow-up 
potentiometer, which is not otherwise required in this use of the VARIPLOTTER. A voltage z may be 
applied to one end of this potentiometer, with the other end grounded. The output voltage at the wiper of 
the potentiometer will then be (z) = f(x). For further details and instructions about the connections of the 




+ 10 



6 



/n [x(0) + a] 



10 



0 ' N v - Y= -/n ( 



x+ a ) 




♦Hr] 



+ (x + a) ~\ -X 
I 




xo 
4 # + x 



ri 9 



a_ 
10 



7045 



4 X - 



-10 



Figure 38. Generation of y = In (x + a) 
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VARIPLOTTER and the curve follower the reader is referred to the basic 1100 VARIPLOTTER manual 
and the Function Generator Model 2.040 manual. 

EXAMPLE Consider a spring-mass-damping system similar to the one discussed earlier. Let us assume 
that the damping coefficient C is now a function of the displacement x. The equation des- 
cribing the system then becomes : 

mx + f (x) x + kx =0 (5 6 ) 

where C - f(x) is given in form of a graph. (Figure 39.) 

LFT/SECj 



~X[ft] 



+ X fjft] 



Figure 39. Graph of Damping Coefficient C 

In the investigation of this system one would prepare the graph for the VARIPLOTTER and curve 
follower as outlined above, and the computer circuit would have the following form (scale factors have 
been neglected). See Figure 40. 

Another device, BIVAR (bi-varient function generator), may be used to generate functions of two 
independent variables. 



o 



- x 




-f(X)X 



/ 



/ 



O 



+ x 



VARIPLOTTER ARM DRIVE 



f(X) 



Figure 40. Function Generation with VARIPLOTTER and Curve Follower 
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TRANSFER FUNCTIONS 

The engineer often writes differential equations describing a system in a certain operational form. This 
engineering procedure uses "Transfer Functions" which are usually defined in terms of "Time Constants" 
These terms may be illustrated by the following simple example: 

Consider a capacitor connected in series with a resistor, switch and battery (Figure 41). At time t after 
the switch is closed, the current flowing in the circuit will be i and conditions are described by tne 
equation : 

E - Ri - i-f i dt = 0 (57) 



which may be readily solved by several methods. It is customary to use an operational notation in such 
equations, and one may write 

[ ] 



/ [ ] dt = 



and 



d[ ] 
dt 



= p[ ] 



(58a) 
(58b) 



Equation 57 may be written in an operational form 



E - Ri - i- i =0 
C P 



(59a) 




T 



TRANSFER FUNCTION 



(WHERE T = RC) 



Figure 41. Simple Circuit and its Transfer Function 
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Whence, by treating the operator p as though it were an algebraic quantity: 



i = (E) (— SL_ ) (59b) 
1 + RCp 



the output voltage V 0 = gr 



P 



so that f. - un ) (60) 

The transfer function of the circuit may be defined as the ratio of the output quantity to the input quantity . 
In this simple case, where the quantities are expressed in the same units, the transfer function G(p) is 
given by 

G(p) - - — — (61a) 
E 1 + RCp 

The quantity RC has dimensions of time and is termed the time constant r (The operator p has dimen- 
sions — - — ) 
time 

One may write in this case : 

C(p) - — L_ (61b) 
1 + rp 

and this function is sufficient to describe the voltage characteristics of the circuit. 

If the transfer functions of a system are known or obtainable, one may readily simulate the system on an 
electronic analog computer. 

The computer circuit for the above example is obtained as follows : 

From equation: V a + ?p V 0 = E (62) 
Solving for the highest derivative : 

p V 0 = - - ^ (63) 
t r 

Assuming that the highest derivative exists at the input to an integrator we obtain the computer circuit 
shown in Figure 42. 

Occasionally the engineer will be confronted by a transfer function representing a system, which is des- 
cribed by a polynomial in p, divided by a polynomial in p, such as: 

E± = K lP 2 + K 2 p + K 3 . 
E { " p* + K 4 p> + K sP + K 6 

The following approach is recommended for the setup of this problem: 
As a first step expand the equation into differential equation form, giving: 

p 3 E 0 + K A p>E 0 + K 5 pE 0 + K 6 E 0 = K^E t + K 2 pE t + K 3 Ei (65) 

One now proceeds by dividing the equation by p"" 1 , where n equals the highest power of p. In this case 
the divisor will be p a . 

pE 0 + K 4 E 0 + K s Eo - + K 6J B 0 4" = K * E i + K > E i - + K > E i \ (66) 
p p 2 p P 
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- VO 



Figure 42. Circuit Representing the Transfer Function G(p) = + r 
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"•"^ (K 3 Ei - K 6 EQ)+^(K 2 Ei- K 5 Eo) 
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£ -K fi E 



6^0 



+ K 3 E I 




+ 7(K 6 Eo-K 3 Ej) 



5 E 0~ K 2 E i 




+ K 2 E- ( 



K2 



Figure 43. Computer Circuit Representing the Transfer Function 



q/ \ _ KiP 2 + K 2 p + K 3 
{PJ p 3 + K A p> + K sP + 



Solving the equation for the highest derivative and grouping the terms according to powers of p gives : 

P E 0 = X^i - K 4 E 0 + -i_ ^,£ f - K,£ 0 ; + ^^jE; - K 5 E 0 ) (67) 

p a p 

If one remembers that the notation ^"represents an integration and if it is assumed that Ei and E 0 are 
known, the equation can then be programmed yielding the computer circuit shown in Figure 43. 
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III. MISCELLANEOUS ADVANCED TECHNIQUES 
REPRESENTATION OF DISCONTINUITIES 

A variety of problems require the simulation of non-linear effects such as displacement, limiting, hys - 
teresis, gear backlash, etc. Most of these functions may be reasonably well represented on the analog 
computer using the signal voltage comparator and/or diodes in conjunction with operational amplifiers. A 
few basic computer circuits commonly used to represent non-linear phenomena are given below. 

1. Simple Limiting 

Desired response: — Z x = X = l 2 (See Figure 44.) 

The amplifier output will follow the input voltage until one of the diodes conducts. This will introduce a 
low impedance into the feedback path. The feedback impedance will not be zero because of the 
attenuator impedance — — - R a tt, and the forward resistance of the diode. Therefore a small gain will 
exist which will give a slight rise to the limited portion of the output voltage. 




Figure 44. Simple Limiting Circuit 

2. Bang- Bang Circuit 

Desired response: ^ limited = ~ ^i f° r ~°°<X<0 

^limited = + h for 0 < Z < °o (See Figure 45.) 

Note that the feedback resistor was omitted in this circuit. The output voltage of the amplifier will there- 
fore very rapidly swing in one direction or the other depending upon the sign of the input — [x] until one 
of the diodes conducts. 

3. Zero Limiting (Half wave rectification) 

Desired response: ^limited ~ 0 for — oo < X <0 

limited = X for 0<X<oo (See Figure 46.) 

The diode at the output of the amplifier will introduce a high impedance when the amplifier attempts to go 
negative. The diode around the amplifier limits the amplifier in the negative direction for positive inputs. 
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-DO 



* V 



+ [x LIMITED ] 



<> + 
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+ [x] 





Figure 45. Bang'- Bang Circuit 



X LIMITED 




X LIMITED 



+ W 



Figure 46. Zero Limiting Circuit 



4. An Absolute Value Circuit 

DesiFed response: X = j X | for all X (See Figure 47.) 

5. Bang- Bang Circuit Using the Signal Voltage Comparator . (See Figure 48.) 

A more extensive tabulation of circuits representing common non-linearities could be easily developed. 
However, with the few basic circuits shown above the reader will be able to develop other circuits suit- 
able to the requirements of a particular problem. 




+ [x] 

Figure 47. Absolute Value Circuit 
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10 



IN 



IN. 



6^ 

V' io 




-10 



- 10 



Figure 48. Bang- Bang Circuit Using the Signal Voltage Comparator 



INTRODUCTION TO THE SOLUTION OF PARTIAL DIFFERENTIAL EQUATIONS 

Partial differential equations can be successfully solved on the analog computer. However, a mathematical 
manipulation of the equation is necessary for the analog computer will only integrate with respect to one 
variable, i.e. time. 

In solving partial differential equations some method must be employed which will transform the partial 
differential equation into a set of ordinary differential equations. By using finite difference approxima- 
tions for the derivatives with respect to all but one of the independent variables a set of finite-difference 
differential equations are written. 

For introductory purposes let us assume we are to investigate the equation for the natural vibration of an 
elastic string : 

d 2 y p d 2 y 
dx 2 ~ T dt 2 



(68) 



where p is the mass per unit length and T is a constant tension. For such an equation the values of y 
and should be known as a function of x for t = 0. These constitute the initial conditions of the prob- 
lem. The boundary conditions must be known at both ends for all values of time t. 

In order to reduce the partial differential equation to a set of ordinary differential equations, let us ex- 
press the space derivatives in terms of finite differences. We select three points as shown in Figure 49. 
spaced Ax apart. We can then obtain the approximate values of the first derivatives at the midpoints of 
the two Ax intervals : 

1 



dx 



Ax 
2 



A* 



(jx - Jx - Ax) 



(69) 



dy 
dx 



x + 



(y* + Ax - y x ) 



(70) 



From equations 69 and 70 it is apparent that the accuracy of the approximation depends upon the size of 
the chosen interval Ax. Dividing the length I into segments Ax effectively means that we are solving the 
problem at discrete points, rather than continuously. 



Let us now consider the second derivative for point x. 
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It if 



d 2 y 
dx 2 



_L t§XL 
Ax \ dx 



x + 



dx 



Ax , 

~2~ 



(71) 



dx 2 



Ax 2 



(y% + Ax - 2 Yx + J% - Ax) 



Thus for position x we can write the equation for the natural vibration of the elastic string 



dt 2 



Ax 2 



hx + a% - 2 y% + j% - Ax) 



Equation 73 can be generalized for any arbitrary point n along the string so that 



d 2 y 
dt 2 



hn + 1 - 2y„ + y n _ i ) 



n pAx 2 

We shall now divide the string into five segments of length Ax = l/ s as shown in Figure 50. 
The boundary conditions are : y 0 = 0 y 5 = 0 

From equation (74) we obtain the following finite difference equations: 

d*y x 25 T . s 

— — = (y 2 - 2yJ 

d* 2 p/ 2 

J 2 y 2 25 r , _ v 

* ■ ^ <* - ^ + J") 



d 2 y 3 25 r 
dt 2 pi 2 



(y 4 - 2y 3 + y 2 ) 



(72) 



(73) 



(74) 



(75) 



(76) 



(77) 
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Figure 50. 



«I ( - 2y . + y .) (78) 

Thus, by dividing the string into n = 5 segments we have changed the partial differential equation to 
(n — 1 ) = 4 second order linear differential equations which can easily be programmed on the computer. 
The unsealed computer circuit is shown in Figure 51. 

The above method for solving partial differential equations on the analog computer is the one most fre- 
quently used. Solutions are usually restricted to simple cases in two or possibly three dimensions. More 
information regarding computing methods for solving partial differential equations can be found in papers 
listed in the bibliography. 

INTRODUCTION TO THE SOLUTION OF ALGEBRAIC EQUATIONS 

Although the primary purpose of the analog computer is its use for the solution of differential equations, 
algebraic relations may have to be handled as well. One might wish to solve a set of simultaneous al- 
gebraic equations, and on the other hand, "algebraic loops" (circuits without integrators) can also ap- 
pear in the computer circuit describing the dynamic behavior of a system. The necessary considerations 
of such algebraic relations are briefly described in the following. 

Consider the set of algebraic equations: 

«n*i + a 12 x 2 = b x (79) 

a 21 x t + a 22 x 2 - b 2 (80) 

where a xl and a 22 are both positive coefficients. If we solve these equations using a similar approach as 
for differential equations, we obtain: 

*i = - °^x 2 (81) 



a., a 



ii 



x 2 = — — — x, (82) 



a„ a 



22 
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Figure 51. 



The computer circuit solving these equations is shown in Figure 52. The circuit might have only two amp 
lifiers in the loop. A condition of positive rather than negative feedback then exists and the circuit will 
be electronically unstable unless the total loop gain K is less than unity. 

K = (a " } (a " } < 1 (83) 
(a lx ) (a 22 ) 



62 



TR- 10 OPERATOR'S MANUAL 




INVERTER IS PRESENT 
IF (a 2| ) IS NEGATIVE. 



INVERTER IS PRESENT 
IF(a |2 ) IS NEGATIVE. 



Figure 52. 



This requirement exists for all algebraic loops consisting of an even number of amplifiers. If it is not 
satisfied, an overloaded condition will quickly arise and no solution of the programmed problem will be 
obtained. Loops consisting of an odd number of summing amplifiers appear to be free of any difficulty of 
this nature. But unfortunately this is not true, for when no integrator is present in a loop with an odd 
number of amplifiers, phase shift unavoidably present at high frequencies within each amplifier will, for a 
given minimum loop gain, produce unwanted oscillations preventing stable conditions from occurring with- 
in the circuit. Thus, no matter the number of amplifiers, in a loop without an integrator one must be very 
careful that the loop gain is small if any problem solution is to obtained. 

For these reasons a naive approach to the solution of algebraic equations is not guaranteed of success. 
A better approach for the setup of algebraic equations is to transform the set of algebraic equations into 
a set of linear differential equations having, as a steady state solution, the solution of the original alge- 
braic equations. Again, considering the above example, we form the set of differential equations: 

«u*i + a 12 x 2 -£>,=- -jj- (84) 



where = 0 in steady state. 
dt 



a 21 x 1 + a 22 x 2 — h 2 = — (85) 

dt 



The computer circuit solving the equations is shown in Figure 53. This technique of solving simultaneous 
algebraic equations, which introduce integrators into the circuit is more acceptable because the integrat- 
ors act as low pass filters which attenuate high frequency oscillations. A stable solution is assured if the 
inverter is present. However, the system may not reach steady state, if the product a 21 a 12 is positive and 
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NVERTER IS PRESENT 
IF (a 21 ) IS NEGATIVE 

INVERTER IS PRESENT 
IF (a |2 ) IS NEGATIVE 



Figure 53. 



no inverter is required. Thus, once more whether or not a solution is obtained depends on the particular 
problem and its programming. Several techniques for the guaranteed solution of algebraic equations with 
the analog computer are known and may be found in papers listed in the bibliography. One such technique, 
that of steepest ascents is briefly described in the following. 

METHOD OF STEEPEST ASCENTS 



This method can be employed when it is desired to find the values of x^ ( i = 1, 2, 
lowing relations are satisfied. 



i = 1 



m) so that the f ol- 
(86) 



where /j are given real valued functions of the x{ and the b{ are constants, and 



dfi f i = 1 .... n 
dxj \ j = 1 .... m 



are continuous. The computer solution is mechanized in such a fashion that the residual errors , de 
fined as fi — b{ (i - 1 .... n), are minimized. In using a least square formulation for the error we 
define n 



X € i 2 
i = 1 



(87) 



and 



whe 



dt 



ds n 

— = 2 X u 
dt i = i 



m 

2 



d e • dx • 
1 J 



det 
dt 



i = 1 



(88) 



(89) 



/=1 dxj dt 
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From equations 88 and 89 one obtains : 



ds " ™ d €i dxj 

dt i = I 7=1 dxj dt 



(90) 



Interchanging summations, equation 90 becomes 



ds m dxj n dei 

— = 2 X 2 (91) 

dt j = 1 dt i = 1 (9x; 



The computer is now mechanized such that : 



dxj dei . . 

— = - a X ££ (92) 

dt dxy 

where a is a positive constant which determines the rate with which the solution converges. From equa- 
tions 91 and 92 one obtains: 



ds _ 2_ ^ / dxj\: 
dt a j = 1 \dt / 



(93) 



Note that the right side of equation 93 will always be negative, i.e. the computer continually tries to de- 
crease s until = 0. The values of Xj change until the minimum value of s is found. Where m = n, 

dt. ' . . . . 

€j- = 0 and a unique solution for the Xj is obtained. If m < n, all n equations cannot be satisfied and 

then a solution is obtained for which s is a minimum. If m < n more equations or constraints are necess- 
ary to determine a solution, for otherwise there are an infinite number of solutions. 

A simple example will illustrate the use of the method of steepest ascents. 

Assume a set of simultaneous algebraic equations of the form 

aa^j + a 12 x 2 = (94) 

(95) 

Rewriting the equations for errors, yields : 

a ll x 1 + a 12 x 2 — = (96) 



With equation 92 we obtain 



a 2X x x + a 22 x 2 — b 2 = e 2 (97) 



= - a Ua^ + e 2 a 21 ] (98) 

dt 



= _ a [fiO u + € 2 a 22 ] (99) 

dt 



Equations 96 through 99 can now be programmed on the computer. The resulting circuit is shown in 
Figure 54. 
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APPENDIX I. SYMBOLS FOR COMPUTER CIRCUITS 



POTENTIOMETERS 



a) GROUNDED 



I ) 0.675 



b) UNGROUNDED 



2 ) 0.500 



OPERATIONAL AMPLIFIERS 



a) INTEGRATOR 



b) SUMMER 



c) HIGH GAIN 




SPECIAL FEEDBACK 
CONNECTIONS 



X 2 DIODE FUNCTION GENERATOR 
TYPE 16.101 



+ X IN 



16.101 



2— I ■ 



10 



MULTIPLIER TYPE 7.045 



+ X I- 



7.045 



+ Y - 



XY 
10 



DUAL FUNCTION SWITCH TYPE 2.127 



up 



dn 



SIGNAL VOLTAGE COMPARATOR TYPE 6.143 



IN, 



+ » • 



IN, 



6.143 



- * • 



* * • A ? 
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Differential Equations 
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"Hold" 

Initial Condition 
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INDEX (continued) 

Subject p age 

Monitoring Panel $j 

Mono tonic Functions 29 

Multiplication by a Constant 8-10 

Multiplication of Variables 75 

Multiplier 75 

Non- Linear Problem 47 

Null Pot 6, 32 

Nulling 5, 31 

"Operate*' 31 

Operational Amplifier 6 

Output Devices 2 

Overload Alarm ', 25, 31 

Partial Differential Equations 59 

Plotters 2 

Pot Set Bus 5, 31 

Potentiometer 5 

Power Switch 31 

Programming 35 

Push -Button , 3 

*'Quarter- Square" X 15 

Recorders x 2 

Rectification (Absolute Value) 58 

Reference Voltage 31 

Repetitive Operation 33 

"Reset" 31 

Scaling 36 

Sign Changing .. m 8 

Signal Voltage Comparator ; 28 

Simultaneous Differential Equations 43 

Square Root 18 

Squaring 17-19 

Static Check 38 

Steepest Ascents 64 

Summation, Algebraic * 12 

Symbols 67 

System Equations 35 

Time Dependence, Explicit 39 

Time Scaling 39, 42 

Transfer Functions 54 

Variable Diode Function Generator 20 

Voltmeter 31 

X 2 Diode Function Generator 17-18 
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